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^ ' Abstract 

I We give a short review of our recent analysis [IJ of the deep inelas- 

tic scattering data (provided by BCDMS, SLAC, NMC) on F2 struc- 
• ture function in the non-singlet approximation with up to next-to-next- 

to-leading-order accuracy and analytic and frozen modifications of the 
strong coupling constant featuring no unphysical singularity (the Lan- 
dau pole). Improvement of agreement between theory and experiment, 
with respect to the case of the standard perturbative definition of Og 
considered recently in [2], was observed. 

> 

in 

2 : 1 Introduction 

' At present an accuracy of the data on the deep inelastic scattering 

(DIS) structure functions (SFs) makes it possible to study separately 
the Q^-dependence of logarithmic QCD-inspired corrections and those of 
power-like (non-perturbative) nature (see for instance [3] and references 
^ . therein). A part of the power corrections may be responsible for the 

^ I deviation of the "physical" strong coupling constant from the standard 

MS one. Similar effect was observed recently ^ in the region of small 
X. In a sense, the study done in [1] is the extension of the analysis [Ij 
to the region of the intermediate and large x values. 

In the present paper we show the results of our recent analyse [T] 
of of DIS SF F2(x,Q2) with SLAC, NMC and BCDMS experimental 
data involved [6]-[B] at NNLO of massless perturbative QCD. This has 
become possible thanks to the results on both the a^{Q'^) corrections 
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to the splitting functions (the anomalous dimensions of Wilson opera- 
tors) [9] and the corresponding expressions of the complete three-loop 
coefficient functions for the structure functions F2 and Fl [IO]. 

As in our previous paper the function F2{x,Q'^) is represented as a 
sum of the leading twist F2^^^ {x,Q'^) and the twist four terms: 

F,(.,g2) = Ff«^^(.,g2)(i + Mj)^ . (1) 

As is known there are at least two ways to perform QCD analysis over 
DIS data: the first one (see e.g. [l3l|T3]) deals with Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) integro-differential equations [T^ and 
let the data be examined directly, whereas the second one involves the 
SF moments and permits performing an analysis in analytic form as 
opposed to the former option. In this work we take on the way in- 
between these two latter, i.e. analysis is carried out over the moments 
of SF F2{x,Q'^) defined as follows 

M^{Q^)= f\'^-'^F^{x,Q^)dx (A; = pQCD, twist2 . . .) (2) 
Jo 

and then reconstruct SF for each by using Jacobi polynomial ex- 
pansion method [16l[T7] (for further details see [5l[18]). The theoretical 
input can be found in the papers [lH [19] . 

2 Infrared modifications of the strong 
coupling constant 

Here, we investigate the potential of modifying the strong-coupling con- 
stant in the infrared region with the purpose of illuminating the problem 
related to the Landau singularity in QCD. 

Specifically, we consider two modifications, which effectively increase 
the argument of the strong-coupling constant at small values, in 
accordance with [20j. 

In the first case, which is more phenomenological, we introduce 
freezing of the strong-coupling constant by changing its argument as 

^ + M^, where Mp is the rho-meson mass [21j (see also [22l [23] 
and references therein). Thus, in the formulas of the previous section 
the following replacement is to be done: 

«i'^(Q')^aP''(Q') = «F(Q' + M2) (i = 0,l,2), (3) 

^For the odd n values, the corresponding coefficients can be obtained by using the analytic 
continuation [TTl [T2] . 
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where the symbol i + 1 marks the ith order of pertmbation theory. 

A second possibihty is based on the idea proposed by Shirkov and 
Solovtsov [231 [25] (see also recent reviews [26j and references therein) 
regarding the analyticity of the strong coupling constant in the complex 
Q^-plane in the form of the Kallen-Lehmann spectral representation. 
This approach leads effectively to additional power Q^-dependence for 
the DIS structure functions (see Eq. ()A2p ). 

2.1 Naive case 

This modification is in a sense quite similar to the freezing procedure 
given above ([3|), i.e. the one- and two-loop coupling constants af'\Q'^) 
and ai^^(Q^) are to be replaced as follows: 

af\Q') ^ a(°)'-(g2) 

a«(Q^) ^ ai^)'^"(Q^) 

where the ellipsis stands for cut terms which give negligible contribu- 
tions in our analysis. 

At the one-loop level, the expression for the analytic coupling con- 
stant dH) is very simple. However, at higher-loop levels, it has a rather 
cumbersome structure (see a recent paper [27j and discussions therein). 
Therefore, it seems to be simpler to use some proper approximations so 
as to be able to carry out a numeric analysis. 

Considering the study ^25j, in the NLO case the difference between 
analytic and standard coupling constants can be represented in the form 
given in ([5]) , which is similar to the LO one with the additional coefficient 
equal to 1/2. Note that numerically this NLO term is quite analogous 
to the LO one, since A(o) <S ^(i) (see, for example, |28j). 

Following the logic expounded in [23], where it was shown that at 
the NNLO level the effective LO Aj^^^ can approximately be taken to 
be 

Al^l = {2n')^A^^A, A^A(2) (6) 



1 ^(0) 
/3o - AJ 

1 



(0) 



2/3o Q2 _ A2 



+ 



(4) 

(5) 



(1) 
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we can apply a simple analytic form in the NNLO as follows H 

a,{Q^)^aT{Q') = a,{Q')-^^^^^+..., as^af\ (7) 

Thus, we propose to use this last expression in the NNLO approxima- 
tion. 

In a sense, the replacement quoted in ([H), ([5]) and d?]) is a naive way 
of doing "analytization" ; we apply the latter procedure to the coupling 
constant itself without considering its actual function. Nonetheless, this 
procedure has already been successfully applied in \12\ U] for analyzing 
the DIS structure functions at small x values. H 

With this motivation we would like to investigate its effect in the 
present study as well. Here it will be referred to as a procedure of the 
naive "analytization". 



2.2 Transition to the canonical form 

To accomplish the procedure of "analytization" more accurately 0, it is 
convenient for the moment M^^{Q^) to be represented in the following 
form (see J19j and discussions therein) 



2/3o 



M-'{Q'') = MriQ'o)-yriQ^Qi)\ («) 

where the new "moment" fi^^{Q'^) starts with as{Q'^) (see [1]). 

The procedure dH) has already been used in [30l[31], where the Grun- 
berg's effective method [32j has been incorporated into the analyses of 
DIS structure functions. 

Now, the "moment" ^^^{Q"^) has the form close to that obtained for 
the sum rule, because it begins with the first power of as(Q^). Conse- 
quently, the form gets closer to that in the difference between the QCD 



^To have the poles in ([5]) and ([7]) exactly cancelled by those in the perturbative expansions 
of QCD coupling (see, for example, j2j), we keep — A"^-^ {i — 0, 1, 2) in the denominators 
of the additional terms (i.e. in the last terms of ([5]) and Q) and, therefore, above LO we 
will have additional terms coincinding with those in the LO case (with the corresponding 
replacements Aq Ai and Aq A) multiplied by the additional factor l/(2z). 

■^The results obtained in [4] show that naive analytization as well as freezing of the 
coupling constant ([3]) leads to the strong improvement of the agreement with experimental 
data for the structure function F2 and its slope at small x values in the double asymptotic 
scaling regime (see and references therein). Moreover, these results are very similar 
to the ones [55] obtained in the framework of the infrared renormalon model [3 for HT 
corrections. 

^We will call this case an ordinary analytic perturbation theory (APT). 
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sum rule and its Parton Model value. Following [23j, the analytical 
version of Eq. ([8]) has the following form [1] 



Af\Q') + A?\Q^)B'}^l{n) 



(9) 



i{3)^Ai2xn{l) 



{3)^r)2^R{2) 
NS 



n] 



where Am is the "analytized" m-th power of i-loop QCD coupling 
and the coeffcients Bj!^g{n) and z^g{n) can be found in |T]. 

Thus, in the APT case the procedure features more complicated 
functions rather than just the powers of some coupling constant. In the 
one-loop case the Euclidean functions of analytic perturbation theory 
(APT) A'j}\Q^) are found to be [25] 



1 

Wo 
1 



1 



L 



(0) 

1 



A 



(1) 



1 dA 



(1) 



1 



a^'Hq') 
aS\q') 

aS\q') 

where L(o) = ln(Q^/(A(o))^), while beyond the LO approximation 
[23] the transform from standard perturbation theory to the APT is 
slightly modified to assume the following form: 



A2 ^3 



(10) 



{4\Q 



A 



A 



A 



1 



4) 



2i/3o g2 _ A2 



1,2) 



, (11) 



2Z/32 (Q2_A2p2 
1 A?.)Q^(Q^ + A?,) 



2^/3,3 



where the NLO and NNLO results (for i 
approximations. 



2(Q2 
1 and 2) are only some 
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To clear up with the meaning of all these formulas, we would like to 
note one more time that in the expressions for Am{Q^) the lower sub- 
script stands for the power of the coupling constant, while the upper one 
is related with the order of the approximation considered. Therefore, if 
the a^-expansion of some variable starts with the first power, as is the 
case at hand, then the power of the last term of the expansion coincides 
with the order of the approximation, i.e. for the last term upper and 
lower subscripts coincide. 



2.3 Fractional analytic perturbation theory 

Recently, in a series of papers [33j) the analytic continuation (jlOp has 
been extended to the noninteger powers of LO L^oJ. In the one-loop 
case, this so-called fractional APT (FAPT) gives: 



Lii_^(e-^(o)) 



(12) 



^(0) ^(0) r(i.) 

where 

Li.(-) = E ^ (13) 

m=l 

is actually the polylogarithm function. H 
It is clearly seen that by virtue of 

Li„„) = ^,LL„(^)=(4)"^ = (4)"_1_, (14) 

all the equations quoted in (llOp can be shown to be well reproduced. 
Since the product of F-functions satisfy 



sin(7rz/) 

it is easy to obtain the following representation (see Appendix A in [Ij): 

1 

^h-u{z) _ z sin(7rt/) f dg ,^-^fl\ 

~T(^ - ^ J 1^^'' [~d ' ^ ^ 

that holds for Ke(l — i^) > and Ke(z^) > (in the case at hand these 
boil down to < < 1). 



^In |33l it was called the Lerch transcendent function. 
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Then, following the previous section let's recast Eq.p2p in the fol- 
lowing fashion: 

Lii_ . (A^o) /Q^ 

which reproduces the equations in (llOh with v = 1,2 and 3, in order. 

Above LO, using quite the same arguments as in the previous sub- 
section we have by analogy (i = 1,2): 



(a<»)(Q=))''^4"«?=) = - ^ (16) 



which in turn reproduces a set of equations quoted in (|lip for v = 1,2 
and 3, respectively. 

Note that the Mellin moments used in this case, contain 
I, = 7j').(n)/(2/3o) + i with i = in the LO, i = 0, 1 — NLO, and i = 
0,1,2 — NNLO approximations. The argument of the polylogarithm 
function A(j)/Q^ found in Eqs. ([I6]) and ([IT]) can be expressed through 
the strong coupling constant (see |lj). Therefore, the power v lies within 
the range < < 4, because < -/Pg{n)/{2f5o) < 2 for 2 < n < 10 
used in the analyses. The integral representation given in Eq. (jl5|) is 
correct only for < < 1, hence the need to extend it to higher values 
of u. Omitting details of this latter extension [fl we have (0 < (5 < 1) 



T{N + 6) \ dzj T{N + 6)r{l - 5) J 1 - \t 

(18) 

To cover an entire range < i/ < 4, we should consider N = 0, 1, 2 
and 3. For = 0, Eq. (jlSp with the corresponding replacement i' ^ 6 
can be used. For > 0, it is straightforward to obtain 

^d\ z _ z /d\'^z _z{l + z(,) 



dzJl-zS, (1-2^2' \ dzJ l-Z^ (1-Z^)3' 

dY z _ zjl+Az^ + z^e) 
^dzj 1-zC ~ (l-zO^ ^ ' 

and make use of the following formulas valid for < 5 < 1: 



Li-5(^) _ z }d^\n-'{l/0 



T{l + 6) T{1 + 6)ril - 6) J (l-z^) 



(20) 







^This is done in Appendix A of [1], where iV = 1, 2, 3. Similar study can be found in the 
recent paper [34] , 
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r(2 + 6) 

U_2-5{z) 

T{3 + 6) 



T{2 + 6)r{l-6) 



r(3 + ,5)r(i-5) 



(1 - z^)^ 
(1 - z^^ 



1 + zC 



(21) 



l + 4z^ + zY(?2) 



i.e. Eqs. ([20]) . (pTI) and (f22]) can be used within the ranges I < < 2, 
2 < 1/ < 3 and 3 < < 4, respectively. 



3 A fitting procedure 

A numeric procedure of fitting the data is described in the previous pa- 
pers [21 [5]. Here we just recall some aspects of the so-called polynomial 
expansion method. The latter was first proposed in [35] and further de- 
veloped in [36j. In these papers the method was based on the Bernstein 
polynomials and subsequently used to analyze data at NLO [37^ [TT] and 
NNLO level \38\ I39j . The Jacobi polynomials for that purpose were 
first proposed and then subsequently developed in [Ml [IT] and used 
in ^301 El [391 EH 

With the QCD expressions for the Mellin moments M^(Q^) analyt- 
ically calculated according to the formula in ([2]), the SF F2(x,Q^) is 
reconstructed by using the Jacobi polynomial expansion method: 

^ max ^ 

n=0 j=0 

where 0"''' are the Jacobi polynomials and a, b are the parameters fitted. 
A condition put on the latter is the requirement of the error minimiza- 
tion while reconstructing the structure functions. 

Since a twist expansion starts to be applicable only above ~ 1 
GeV^ the cut > 1 GeV^ on the data is applied throughout. 

MINUIT program [41] is used to minimize two variables 



tear 



Xslope 



j^exp _ J) 



tear 



where D = dlnF2/dlnln(5^. The quality of the fits is characterized by 
X^/DOF for the structure function F2. Analysis is also performed for 
the SF slope D that serves the purpose of checking the properties of fits 
(for more details see [2]). 

We use free normalizations of the data for different experiments. For 
a reference set, the most stable deuterium BCDMS data at the value of 
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the beam initial energy Eq = 200 GeV is used. With the other datasets 
taken to be a reference one the variation in the results is still negligible. 
In the case of the fixed normalization for each and all datasets the fits 
tend to yield a little bit worse x^; just as in the previous studies. 

4 Results 

Since there is no gluons in the nonsinglet approximation the analysis 
is essentially easier to conduct, with the cut imposed on the Bjorken 
variable x > 0.25 where gluon density is believed to be negligible. 

Here we conduct separate and combined analyses of SLAG, BCDMS 
and NMC datasets for hydrogen and deuterium targets. The cut on x is 
imposed in a combination with those placed on the y variable as follows: 



y 


> 0.14 


for 


0.3 < a; < 0.4 


y 


> 0.16 


for 


0.4 < X < 0.5 


y 


> 0.23 


for 


0.5 < X < 0.6 


y 


> 0.24 


for 


0.6 < X < 0.7 


y 


> 0.25 


for 


0.7 < X < 0.8 



which are meant to cut out those points with large systematic errors. 
Thus, upon imposing the cuts a complete dataset consists of 327 points 
in the case of hydrogen target and 288 — deuterium one. The starting 
point of the QCD evolution is taken to be Qq = 90 GeV^. This Qq value 
is close to the average values of spanning the corresponding data. 
Prom earlier studies [21 [5] it follows that it is enough to take the maximal 
value of the number of moments to be accounted for N^ax = 8 [17j; also 
note that the cut 0.25 < x < 0.8 is imposed everywhere. 

To reduce the number of parameters, we perform two groups of 
fits. The first one is dealt within the variable-fiavor- number-scheme 
(VFNS) [2] and the H2 and D2 experimental datasets analyzed simul- 
taneously. The results for the second group are obtained within the 
fixed-fiavor-number-scheme (FFNS) with an active number of flavors 
n/ = 4 and the H2 and D2 experimental datasets considered separately. 

4.1 VFNS case 

All the results obtained within our reference VFNS [2] and in the cases 
of "naive", APT, and "frozen" modifications of are gathered in a set 
of tables separately for each order of perturbation theory approximation 
and displayed in Figs. 1-3 separately for all these three cases. 
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Table 1. LO values of the twist-four term hi{x) (with statistic 
errors given) obtained in the analysis of the combined H2 + D2 dataset 
within the VFNS and with various modifications of ag 



X 


Naive cmal}^. as 


APT-inspired as 


Frozen as 


Standard as 


0.275 


-0.121 ± 0.008 


-0.123 ± 0.008 


-0.204 ± 0.011 


-0.271 ±0.012 


0.35 


-0.055 ± 0.007 


-0.055 ± 0.007 


-0.167 ± 0.017 


-0.257 ±0.017 


0.45 


0.119 ± 0.012 


0.119 ± 0.012 


-0.021 ± 0.031 


-0.144 ±0.030 


0.55 


0.422 ± 0.022 


0.422 ± 0.023 


0.211 ± 0.053 


0.051 ±0.049 


0.65 


0.870 ± 0.060 


0.866 ± 0.059 


0.558 ± 0.095 


0.364 ±0.088 


0.75 


1.322 it 0.117 


1.336 ± 0.112 


0.917 ± 0.152 


0.709 ±0.138 


X^/DOF 


0.93 


0.93 


0.91 


0.94 




2.30 


2.35 


2.15 


2.60 




0.1474 


0.1474 


0.1409 


0.1400 


Table 2. 


NLO values of the twist-four term /i4(x) ( 


with statistic 


errors given ) 


obtained in the analysis of the combined H2 


+ JJ2 dataset 


within the VFNS and with various modifications of ag 




X 


Naive analyt. as 


APT-inspired as 


Frozen as 


Standard as 


0.275 


-0.127 ± 0.009 


-0.129 ± 0.007 


-0.183 ± 0.008 


-0.229 ± 0.010 


0.35 


-0.098 ± 0.007 


-0.024 ± 0.009 


-0.149 ± 0.010 


-0.218 ± 0.016 


0.45 


0.014 ± 0.012 


0.187 ± 0.013 


0.010 ± 0.019 


-0.084 ± 0.030 


0.55 


0.172 ± 0.024 


0.506 ± 0.019 


0.215 ± 0.033 


0.098 ± 0.052 


0.65 


0.339 ± 0.057 


0.910 ± 0.045 


0.476 ± 0.065 


0.356 ± 0.093 


0.75 


0.478 ± 0.107 


1.230 ± 0.090 


0.757 ± 0.108 


0.648 ± 0.145 


X^DOF 


0.85 


0.84 


0.97 


1.02 


xlioJDOF 


0.82 


0.78 


0.87 


1.20 




0.1275 


0.1224 


0.1169 


0.1152 



Table 3. NNLO values of the twist-four term h4{x) (with statistic 
errors given) obtained in the analysis of the combined H2 + D2 dataset 
within the VFNS and with various modifications of as 



X 


Naive analyt. Qs 


APT-inspired as 


Frozen as 


Standard as 


0.275 


-0.171 ± 0.006 


-0.196 ± 0.008 


-0.149 ± 0.006 


-0.173 ± 0.017 


0.35 


-0.160 ± 0.008 


-0.152 ± 0.012 


-0.129 ± 0.013 


-0.094 ± 0.020 


0.45 


-0.044 ± 0.018 


0.030 ± 0.022 


-0.007 ± 0.031 


-0.110 ± 0.015 


0.55 


0.085 ± 0.033 


0.269 ± 0.038 


0.116 ± 0.062 


-0.086 ± 0.033 


0.65 


0.221 ± 0.065 


0.551 ± 0.074 


0.218 ± 0.115 


0.085 ± 0.083 


0.75 


0.304 ± 0.100 


0.782 ± 0.116 


0.258 ± 0.169 


0.158 ± 0.105 


X'/DOF 


0.98 


1.02 


0.97 


0.92 




1.22 


1.17 


1.02 


1.83 


as(M|) 


0.1151 


0.1125 


0.1163 


0.1159 



From Tables 1-3 it is seen that in all the cases considered /DOF 
shows good agreement between the experimental data and theoretical 
predictions for the SF Mellin moments. In LO and NLO cases, the an- 
alytic and "frozen" modifications lead to some additional improvement 
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of fits, namely, /DOF in these cases is less than that in the standard 
case. 

Since the quantity xiiope/ inherently linked with the pQCD 
aspects to be observed in the data it in this respect is very informative, 
for it strongly varies from one approximation to the other thus indi- 
cating if there are any effects incompatible with the dependence in 
each x-bin assumed. It is seen that it, much like the x^/DOF quantity, 
demonstrates similar tendency: in all the cases considered, the analytic 
and "frozen" modifications lead to improvement of fits. Certain im- 
provement of the quality of fits is observed at the NLO level. In this 
approximation even the standard as case leads to reasonable agreement 
for the slopes; furthermore, for all the infrared modifications it is seen 
that x^gi^p^l DOF < 1. In the NNLO approximation the numbers for the 
slope in the case with a standard as are not as good but the infrared 
modifications (especially the "frozen" one) lead to better agreement be- 
tween theoretical and experimental results for the former quantity. It is 
difficult to pin down a reason for such a deterioration of this agreement 
(in the scheme with a standard ag) when NNLO corrections are added. 
Because this effect is absent in the FFNS case (see the following sub- 
section), we suppose that the exact equation for the coupling constant, 
used here, is somehow in inconsistency with the NNLO expression for 
the heavy quark thresholds, which is in fact based on certain expan- 
sions of the coupling constant at the threshold crossing points (see, for 
example, [2]). We plan to study this fine effect elsewhere. 

The QCD coupling decreases from LO through NNLO, which is in 
perfect agreement with other studies (see [19] and references therein); it 
is seen that the frozen modification gives the results closest to those for 
the standard ag. In the analytic cases the values of as(M^) are higher 
in the first two orders of perturbation theory, with the maxima in 
the difference Aas(M|) being 0.0123 and 0.0072 for the naive analytic 
and APT cases, respectively, observed at NLO. Thus, the differences 
are substantially greater than the values of the total experimental er- 
ror A*°*'^'as(M|) = 0.0022, which was obtained in [2] by combining 
statistical and systematic errors in quadrature. 

At the NNLO the APT procedure leads to lower as(M|), though 
in the naive case the a<j(M^) value is closer to those obtained for the 
frozen and standard versions. Nevertheless, all the NNLO as{M'^) val- 
ues, except for the APT case, are in good agreement within statistical 
errors, which were found to be A'**"*a^(M|) = 0.0007 m our previous 
studies [2]. The APT-inspired QCD coupling constant is compatible 

^This observation is consistent with earlier studies (see [12j and references therein). 
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with the rest within the total experimental error. 

From Tables 1-3 and, particularly, Figs. 1-3 it is seen that in the 
cases of analytic and "frozen" coupling constants the twist-four correc- 
tions are larger compared to those in the case of a standard perturbative 
coupling constant, thus confirming the results obtained in [23]. For ex- 
ample, at X ~ 0.75 the higher- twist corrections (HTCs) in the standard 
case are about twice as less that those obtained in the cases of analytic 
modifications in all orders considered. In the "frozen" case, the HTC 
values are compatible with those in the standard case within statistical 
errors. H 

The difference in the values of HT parameter for the naive analytic 
and frozen variants becomes moderate at the NLO level unlike the APT 
case. For the infrared modifications the HT terms are large and in the 
frozen and APT cases are compatible with the LO ones. A partial expla- 
nation of this effect can be found in Appendix A. The "analytization" 
of the coupling constant generates additional power-like contributions 
with the opposite sign as compared to the twist-four corrections at large 
X and, thus, increases effectively the values of HTCs. The analysis in 
Appendix A is given at the LO level. Above the latter the correspond- 
ing results, in accordance with Eqs. (jH, ([5]) and ([7]), can be estimated 
by replacing A(o) A(i)/2 in NLO and A(o) A/4 in NNLO in front 
of the term ln(l — x) in ()A2p . 

In the NNLO approximation the situation changes drastically. HT 
terms for all the cases, except for the APT case, are comparable with 
each other. HTCs in the latter case are still higher but they are also 
strongly suppressed compared to HTCs in the LO and NLO cases. To 
some extent, this can be explained by replacing A^q^ — )• A/4 in front of 
ln(l — x) in ()A2p . which decreases the influence of the infrared modifi- 
cations. 

Thus, at the NNLO level the twist-four corrections appear to be 
small for all the cases considered above. 

The values for parameters in the parameterizations of the parton 
distributions for the cases corresponding to different coupling constant 
modifications are given in Table 4 of [Ij. 

^For the results to be discernible. Figs. 1-7 contain statistical errors for the infrared 
modifications only. The magnitude of the errors is on par in the standard case and is not 
shown. 
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4.2 FFNS case 



For comparison let's present the values of PDF parameters obtained 
within the FFNS (nj = 4) in the analyses of the hydrogen and deuterium 
data for the versions of the strong coupling constant discussed above, 
with the addition for the case of FAPT- inspired modification of . Here 
we have no threshold transitions for QCD coupling and PDF Mellin 
moments and, hence, are able to consider the hydrogen and deuterium 
datasets separately. 

The values for parameters in the parameterizations of the parton 
distributions for the cases corresponding to different coupling constant 
modifications are given in Table 5 of [IJ. 

In order to be able to assess the difference among the cases with 
frozen, FAPT and standard versions for the strong coupling constant 
let's present the tables with the HT values obtained within FFNS {nf = 
4) for the hydrogen data (results obtained for the deuterium data can 
be found in [Ij). 

Table 4. LO values of the twist-four term obtained within the FFNS 
(uf = A) and with various modifications of Us . Only statistical errors 
are given. 





FAPT-inspircd 


Frozen 


Standard 


X 


h4{x) 


h4(x) 


h4(x) 


0.275 


-0.221±0.011 


-0.183 ± 0.012 


-0.235±0.012 


0.35 


-0.187±0.014 


-0.160 ± 0.018 


-0.232±0.021 


0.45 


0.002±0.023 


-0.023 ± 0.037 


-0.130±0.038 


0.55 


0.332±0.034 


0.189 ± 0.065 


0.049±0.065 


0.65 


1.001±0.063 


0.610 ± 0.117 


0.455±0.111 


0.75 


2.031±0.131 


1.177 ± 0.207 


1.003±0.182 


X^/DOF 


0.98 


0.94 


0.98 


^slope^ 


1.47 


1.25 


1.58 




0.1394 


0.1387 


0.1376 



Table 5. NLO values of the twist-four term obtained within the 
FFNS (nf = A) and with various modifications of ag.Only statistical 
errors are given. 





FAPT-inspired Ob 


Frozen as 


Standard a^ 


X 


h4{x) 


h4(x) 


h4(x) for H2 


0.275 


-0.205±0.012 


-0.154 ± 0.011 


-0.200±0.012 


0.35 


-0.194±0.017 


-0.147 ± 0.016 


-0.219±0.019 


0.45 


-0.067±0.030 


-0.058 ± 0.037 


-0.169±0.041 


0.55 


0.160±0.044 


0.069 ± 0.067 


-0.078±0.072 


0.65 


0.635±0.071 


0.317 ± 0.118 


0.153±0.123 


0.75 


1.512±0.130 


0.723 ± 0.204 


0.534±0.205 


X^/DOF 


0.91 


0.89 


0.92 


xlo,JDOF 


1.08 


0.92 


1.23 


«4M|) 


0.1220 0.1200 


0.1192 
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Table 6. NNLO values of the twist-four term obtained within the 
FFNS (rif = 4j and with various modihcations of Og. Only statistical 
errors are given. 





FAPT-inspircd as 


Frozen as 


Standard as 


X 




h4{x) 


h4{x) 


0.275 


-0.168±0.010 


-0.119 ± 0.011 


-0.158 ± 0.020 


0.35 


-0.181±0.015 


-0.112 ± 0.010 


-0.166 ± 0.021 


0.45 


-0.130±0.037 


-0.049 ± 0.018 


-0.156 ± 0.036 


0.55 


-0.056±0.065 


0.007 ± 0.031 


-0.157 ± 0.079 


0.65 


0.126±0.111 


0.106 ± 0.058 


-0.061 ± 0.127 


0.75 


0.419±0.178 


0.240 ± 0.116 


0.049 ± 0.209 


X^/DOF 


0.89 


0.87 


0.89 




0.97 


0.68 


0.97 




0.1183 


0.1180 


0.1176 



Just like in the previous subsection it is seen from Tables 4-6 that in 
all the cases at hand, we have good agreement between the experimental 
data and theoretical predictions. Once again, in all the cases considered, 
excluding LO FAPT, the analytic and frozen modifications lead to slight 
improvement of fits, that is /DOF and xliope/ found to be 
smaller in these cases than those for the standard q^, and the QCD 
coupling constant decreases when we move from LO through NNLO. 

Also note that contrary to what was observed in the previous subsec- 
tion here the quantity xliope/ ^^-^ steadily decreases when we proceed 
step by step from LO to NLO and then to NNLO level. Respectively, 
here the smaller values of xiiope/ observed for the "frozen" 
case as well. 

The values of coupling constants are very similar, especially for the 
"frozen" and standard versions. In the analytic cases, the central values 
of as(M^) are little higher, mostly in the first two orders of perturbation 
theory (see also [12]). 

In the NNLO, all the as(M^) values, except for the analytic one, are 
in good agreement within statistical errors. As earlier, the analytic QCD 
coupling is in agreement with the rest within the total experimental 
error. 

Similarly to the previous subsection we note in Figs. 4, 5, 6 appre- 
ciable difference between HTCs obtained in the standard and analytic 
cases in the first two orders of perturbation theory. The difference in 
the cases of a "frozen" and standard Og is not as large; the values of HT 
terms are in agreement within statistical errors (see Fig. 7). 

In the NNLO, the situation is over again changed considerably. In 
all the cases of infrared modifications considered, the HTCs are small. 
Similarly to the previous subsection they are not compatible with zero 
at j; ~ 0.75, while the HT terms in standard QCD are compatible with 
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zero and, at the same time, in agreement with all the cases considered 
within statistical errors. 

Then, at NNLO we can see agreement between standard QCD and 
its infrared modifications for QCD coupling as(M^), as well as for the 
respective HTCs (as a rule) within statistical errors. 

5 Conclusions 

The pattern of separating perturbative QCD and HT corrections may 
be different in different orders of perturbation theory, as well as in some 
resummations based on several first orders, and for certain modifica- 
tions of the strong coupling constant as well. In the present paper, we 
have shown the results of our recent analysis [1] , where we studied the 
consequences of the infrared modifications of the QCD coupling con- 
stant — the so-called "frozen" and "analytized" versions. In the last 
case three different options were considered: 

• a simple modification of the strong coupling constant [2^ without 
rearrangement of a perturbation series; 

• an application of the ordinary analytic perturbation theory (see [23]) 
to the "moments" ix^^{Q'^) given in ([8]); 

• impact of the fractional analytic perturbation theory |33j applied 
directly to the Mellin moments M^^ {Q"^). 

To test all these modifications, the Jacobi polynomial expansion method 
developed in [161 HZ] was used to perform analysis of Q^-evolution of the 
DIS structure function F2 by fitting all existing to date reliable fixed- 
target experimental data that satisfy the cut x > 0.25. To the best of 
our knowledge, the study pLj is the first application of the FAPT results 
to the fits of the DIS structure functions. 

The main conclusions are as follows. In the first two orders of per- 
turbation theory, the coupling constant as(-^D rises for all versions 
of the analytic perturbation theory. This observation is in complete 
agreement with recent studies in [42j . 

An increase in the central values of as(M^) is smaller in the NNLO 
approximation, much like in the case of the "frozen" version of a^. 
Nevertheless, in all the cases the as{M'^) values are in agreement mostly 
within total experimental errors. For the NNLO case, the results are 
as a rule compatible between each other within statistical errors. Also, 
we note that within statistical errors only (in our case A^^°'^as{M'^) = 
0.0007 [2]) all the NNLO results for a,(M|) in the FFNS case (see 
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Table 6) are too compatible with the world average value for the coupling 
constant presented in the review [33] |£|: 

a,(M|) = 0.1184 ±0.0007. 

It is also observed that there is some rise of the twist-four corrections 
in the first two orders of perturbation theory, particularly for all versions 
of the analytic one. This observation is in complete agreement with 
recent studies [23j of the Bjorken sum rules, where it was shown that 
there is a reduction of higher HT terms, starting with the twist-six 
ones. Unfortunately, in [T| we were not able to study the twist-four 
and twist-six corrections simultaneously since their contributions are 
strongly correlated. 

In the NNLO, all the cases of infrared modifications of the QCD 
coupling feature nonzero although rather small twist-four corrections. 
Whereas in the case of the standard QCD approach, the HT terms are 
close to zero at large values of the Bjorken variable x. However, the 
NNLO HTCs for all the cases considered are compatible between each 
other within statistical errors. 

In principle, the main difference between the cases with a standard 
QCD coupling constant and its analytic and "frozen" modifications is 
in the strong suppression of the higher twist corrections in NLO and 
NNLO orders of perturbation theory, respectively. 

What is interesting to look for further in the study is the consid- 
eration of the combined nonsinglet and singlet analyses using the DIS 
experimental data within an entire x region, as well as an application 
of certain resummation-like Grunberg effective charge methods [32] (as 
was done in [31] in the NLO approximation) and the "frozen" [21j and 
analytic [25] versions of the strong coupling constant (see [U [22t [23] for 
recent studies in this direction). The effect of N^LO corrections would 
also seem to be important to account for in the subsequent investiga- 
tions, as well as an extension of the FAPT model results for to the 
VFNS case. 

^It should be mentioned that this analysis was carried out over the data coming from 
the various experiments and in different orders of perturbation theory, i.e. from NLO up to 
N^LO. 

^"There are a lot of "frozen" versions of the strong coupling constant (see, for example, 
the list of references in [22j). 



16 



6 Acknowledgments 



The work was supported by the RFBR grant No.l0-02-01259-a. We are 
grateful to S.V. Mikhailov, D.B. Stamenov and O.V. Teryaev for useful 
discussions. 



7 Appendix A 



Consider the large x asymptotic [44] of SF F2{x,Q'^) (for simplicity we 
restrict ourselves to the LO approximation) 



where 



F2(x,Q2)^(l_^)f(4'''(Q^)), 



6(ai°)(Q2)) =6o_Jin(aW(Q2)), d 



(Al) 
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and bo is some constant which can be obtained from the quark counting 
rules [is] . 

The Q^-dependent part in the r.h.s. of (|Aip can be represented as 
follows: 



(l-x) 



-Jln(4"){Q2)) 



-dln(l-x) 



Performing the "analytization" procedure and using 1^ , we have 

^ -(Jln(l-x) 



-dln{l~x) 



(0). 



1 + 



dln(l — x) 



(0) 



(A2) 



From this expression it follows that the power corrections, which were 
generated by the "analytization" of the coupling constant, have the op- 
posite sign as compared to the twist-four corrections at large x and, 
moreover, demonstrate a different asymptotic behavior. These correc- 
tions (taken with the additional sign "— ") increase like ln(l — x) while 
the twist-four corrections behave like 1/(1 — x) (see ^7] and discussions 
therein). This difference in sign between these corrections leads to the 
larger twist-four corrections in the present analysis than those given 
in [2]. 



^^This analysis is carried out at rather large values used in the paper. For low 
consideration can be found, for example, in [46j . 
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Figure 1: Comparison of the HTC parameter h^lx) obtained in LO, NLO 
and NNLO for hydrogen data (the bars stand for statistical errors) between 
our reference VFNS with a standard perturbative [2] and that with naive 
analytic a^. 
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Figure 2: Comparison of the HTC parameter h4^{x) obtained in LO, NLO and 
NNLO for hydrogen data between a VFNS and that with APT- inspired a^. 
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Figure 3: Comparison of the HTC parameter h^lx) obtained at LO, NLO and 
NNLO for hydrogen data within a VFNS between the cases with a standard 
perturbative and "frozen" a^. 
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Figure 4: Comparison of the HTC parameter h^lx) obtained in LO, NLO and 
NNLO for hydrogen data within a FFNS {nj = 4) between the cases with a 
standard perturbative and naively analytized a^. 



25 



1.5 



0.5 







LO a. 



LO APT a: 



an 
s 

NLO a. 



L Q 



an 



NLO APT a;^ 
NNLO^a. 



NNLO APT a: 



an 



0.3 



0.4 



0.5 



0.6 



0.7 



0.8 



Figure 5: Comparison of the HTC parameter h^lx) obtained in LO, NLO and 
NNLO for hydrogen data within a FFNS {nj = 4) between the cases with a 
standard perturbative and APT- inspired a^. 
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Figure 6: Comparison of the HTC parameter h^lx) obtained at LO, NLO and 
NNLO for hydrogen data within a FFNS {nj = 4) between the cases with a 
standard perturbative and FAPT-inspired a^. 
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Figure 7: Comparison of the HTC parameter h^lx) obtained at LO, NLO and 
NNLO for hydrogen data within a FFNS {nj = 4) between the cases with a 
standard perturbative and "frozen" a^. 
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